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1  Introduction  and  problem  formulation 
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estimate  polarisation  properties  in  dielectric  material  targets.  In  ordei  to  accomplish 
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one  must  have  a  viable  model  for  acoustic  pressure-dependent  material  polarization.  In  Sec¬ 
tion  2,  we  present  such  a  polarization  model  founded  on  accepted  mechanisms.  Our  model  is 
based  on  the  Debye  model  for  orientational  polarization.  We  modify  this  model  to  incorpo¬ 
rate  pressure-dependence  and  allow  for  the  electromagnetic/acoustic  interaction  that  results 
in  the  reflections.  Other  approaches  to  modeling  electromagnetic/acoustic  interaction  are 
outlined  in  [10].  In  Section  3,  we  develop  a  theoretical  framework  for  the  inverse  problem 
which  underlies  the  interrogation  technique.  Finally  in  Section  4,  we  provide  sample  inverse 
problem  calculations  to  support  the  feasibilty  of  the  technique. 

Under  certain  assumptions  (see  also  [4]),  including  material  homogeneity  in  directions  per¬ 
pendicular  to  the  direction  of  electromagnetic  wave  propagation  and  the  use  of  a  polarized 
planar  electromagnetic  wave,  we  may  model  the  interrogation  technique  using  a  one  di¬ 
mensional  form  of  Maxwell’s  equations  with  temporally  and  spatially  varying  coefficients. 
Provided  that  these  assumptions  hold,  Maxwell’s  equations  can  be  written 


dz 


H  =  —D  +  aE+Js 


(1) 

(2) 


where  E  and  H  are  the  electric  and  magnetic  fields,  D  is  the  electric  flux  density,  and 
is  the  source  current  density.  The  macroscopic  polarization  P  is  introduced  via  the  electric 
flux  density  through  the  relation 


D  —  CoCrE  +  P  (3) 

where  represents  the  effects  of  instantaneous  polarization  in  the  dielectric.  We  discuss  the 
form  of  P  in  great  detail  in  Section  2.  By  taking  the  appropriate  derivatives  of  (1)  and  (2) 
and  employing  equation  (3),  we  obtain 


IXQtrE  4-  p-qP  +  jJLoaE  —  E"  =  —fiojg.  (4) 

Here  and  throughout  we  use  E  to  denote  f^E  and  E'  to  denote. 

We  consider  (4)  in  the  domain  0  <  z  <  1  with  air  in  the  interval  [0,  zi]  and  a  dielectric 
material  in  the  region  (0  <  2:1  <  1).  We  enforce  an  absorbing  boundary  condition 

on  the  left  {z  =  0)  and  a  perfectly  conductive  boundary  condition  on  the  right  (z  =  1). 
We  leave  the  initial  conditions  for  the  electric  field  in  general  form,  but  we  assume  that  the 
polarization  and  its  first  time  derivative  are  initially  zero  (this  may  be  done  without  loss  of 
generality).  The  geometry  of  the  problem  is  shown  in  Figure  1. 
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2  Polarization  model 


In  this  section,  we  piesent  and  motivate  our  model  for  pressure-dependent  polarization.  We 
begin  by  discussing  polarization  in  general  and  then  explain  how  it  pertains  to  our  problem. 

2.1  Mechanisms  of  polarization 

Electric  polarization  is  by  definition  the  electric  dipole  moment  per  unit  volume.  The  for¬ 
mation  of  these  electric  dipoles  can  be  caused  by  several  mechanisms  [2],  [3]  which  we  briefly 
summarize  here. 


Electronic  polarization/  Optical  polarization/  Induced  polarization  An  applied  field 
displaces  the  electron  cloud  center  of  an  atom  with  respect  to  its  nucleus.  This  in¬ 
duces  a  dipole  moment.  Electronic  polarization  is  found  in  both  materials  that  possess 
molecules  with  large  dipole  moments  (polar  materials)  and  those  that  do  not  (nonpolar 
materials). 

Atomic  polarization/  Ionic  polarization/  Molecular  polarization  An  applied  elec¬ 
tric  field  may  displace  the  atoms  in  the  molecules,  changing  the  distance  between 
the  atoms,  and  thus  changing  the  dipole  moment.  Atomic  polarization  only  occurs  in 
polar  materials. 
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Orientational  polarization/  Dipole  polarization  Without  an  applied  field,  a  polar  ma¬ 
terial  possesses  permanent  dipole  moments  that  are  randomly  oriented.  When  a  field  is 
applied,  these  dipoles  align  themselves  with  the  field.  Since  orientational  polarization 
is  reliant  upon  the  existence  of  permanent  dipole  moments,  it  only  is  found  in  polar 
materials. 

Interfacial  polarization  The  impurities  and  defects  in  crystal  can  impede  the  flow  of 
charge  created  by  an  applied  field.  The  resulting  charge  accumulation  can  result  in  a 
dipole  moment.  This  type  of  polarization  is  found  only  in  crystals. 


The  multiple  names  for  each  type  of  polarization  can  be  confusing,  especially  when  comparing 
the  research  of  different  contributors.  We  attempt  to  refer  to  each  mechanism  by  the  first 
name  given  above.  We  point  out  that  in  addition  to  this  terminology  both  atomic  and 
electronic  polarization  are  sometimes  referred  to  as  distortional  polarization  [2]. 

In  a  given  material,  polarization  can  be  the  result  of  one  or  more  of  these  four  mechanisms. 
We  are  primarily  interested  in  materials  that  contain  a  high  water-content,  such  as  living 
tissue  or  mud,  so  here  we  focus  on  polar  liquids.  The  polarization  in  this  class  of  liquids 
tends  to  depend  mostly  on  the  orientation  of  permanent  electric  dipoles  in  the  molecules 
(orientational  polarization)  and  the  distortion  of  the  molecules  by  an  applied  electric  field 
(electronic  and  atomic  polarization)  [16].  With  this  in  mind,  we  focus  on  these  polarization 
mechanisms  in  the  remainder  of  our  discussions. 

In  the  presence  of  most  applied  electric  fields,  the  polarization  of  a  high  water-content  liquid 
is  both  distortional  and  orientational.  At  high  (optical)  frequencies  however,  the  electric 
field  oscillates  so  rapidly  that  it  does  not  hold  any  orientation  long  enough  for  the  dipoles 
to  align  with  it.  Thus,  the  orientational  polarization  is  virtually  insignificant  [26].  This 
implies  that  at  sufficiently  high  frequencies,  the  only  contribution  to  the  dielectric  constant 
or  optical  index  of  refraction  is  from  electrical  distortion  [16]. 

Since  the  polarization  of  a  polar  liquid  has  multiple  mechanisms,  we  expect  that  a  complete 
model  must  incorporate  them  all.  Orientational  polarization  is  suggestive  of  a  mechanism 
with  an  exponential  decay  factor,  such  as  the  one  in  the  model  proposed  by  Debye  [2]. 
However,  a  system  rarely  conforms  exactly  to  the  model  described  by  the  Debye  dispersion 
equations  due  to  the  fact  that  the  polarizational  decay  may  not  be  represented  accurately 
by  a  mechanism  with  one  relaxation  time  [17].  On  the  other  hand,  distortional  polarization 
causes  charges  to  behave  somewhat  like  linear  harmonic  oscillators;  thus  it  is  reasonable  to 
model  them  as  such  (the  Lorentz  model  is  an  example).  Neither  of  these  types  of  models  alone 
will  be  sufficient  to  completely  describe  the  polarization  of  a  polar  liquid.  Nonetheless  as  a 
first  step  and  to  illustrate  our  ideas,  we  base  our  model  on  the  Debye  model  for  orientational 
polarization.  Future  modeling  attempts  will  require  systems  of  more  complexity. 
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2.2  The  Debye  model 

The  Debye  model  [4]  can  be  represented  by  the  first  order  ordinary  diflferential  equation 


tP  +  P  =  e.o{e,  -  Coo)E, 


(5) 


or  by 


P{t,z)=  I  9it-s,z)Eis,z)ds 
Jo 

with  kernel 

«(()  =  CX'P  ‘A^'-^oc) 

In  these  equations,  Cj  is  the  static  relative  perrnitivitty  and  is  the  value  of  permittivity 
for  an  extremely  high  infinite)  frequency  field.  In  this  model,  the  value  of  the  relative 
permittivity  e,.  of  (3)  in  the  dielectric  is  given  by  e^;  that  is,  =  1  in  [0,  zj]  and  in 

(zi,l].  The  variable  r  is  the  relaxation  time  of  the  dielectric. 

In  [2],  Anderson  describes  a  potential  double  well  formulation  for  an  atomic  model  that  leads 
to  the  Debye  polarization  model.  In  this  model,  the  dielectric  is  made  up  of  independent 
noninteracting  particles;  eacli  particle  has  two  equilibrium  positions  separated  by  a  barrier 
of  high  potential.  One  considers  a  charged  particle  with  two  equilibrium  positions  A  and 
B  located  a  distance  d  from  each  other.  Between  them  is  a  potential  barrier  W  such  that 
>  ksT  where  /cb  is  the  Boltzmann  constant  and  T  is  the  temperature.  (See  Figure  2.) 
If  there  is  no  electromagnetic  field  present,  one;  assumes  that  the  particle  oscillates  about 
either  equilibrium,  and  on  occassion,  obtains  enough  energy  to  cross  the  potential  barrier 
and  jump  into  the  other  well.  Over  time,  for  constant  temperature,  the  particle  is  near  A  as 
often  as  near  B  and  the  probability  of  finding  the  particle  near  a  given  well  is  |. 
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d  d 

Figure  2:  Potential  double  well  model  with  and  without  an  applied  field 


When  an  electric  field  E  is  applied  in  the  direction  from  A  to  B,  the  potentials  at  each 
equilibrium  are  no  longer  equal,  for  instance  Va  >  Vb,  and 


14  -Vb  =  edE, 


where  e  is  the  charge  of  the  particle.  (See  Figure  2.)  A  result  from  Boltzmann  statistics 
implies  that  the  probability  that  a  particle  has  potential  V  is  proportional  to  exp(-^), 
so  that  now  it  is  more  likely  to  find  the  particle  near  equilibrium  B.  As  before,  a  particle 
can  jump  from  one  equilibrium  to  the  other  if  it  acquires  enough  energy.  For  a  potential 
barrier  W,  the  probability  that  a  particle  can  cross  this  barrier  in  the  direction  from  B  to 
A  is  proportional  to  exp(— with  proportionality  constant  the  assumed  frequency 
of  oscillation  due  to  thermal  agitation  of  the  particle  about  the  equilibrium.  Likewise,  the 
probability  that  the  particle  can  cross  the  barrier  in  the  direction  from  A  to  B  is  given  by 
^  exp(-fi^^^).  Using  these  probabilities  and  the  fact  that  the  total  number  N  —  Na  +  Nb 
of  particles  is  constant,  one  can  derive  [2]  (see  also  page  387  of  [15])  a  linear  first  order 
differential  equation  to  describe  the  difference  A^b(0  -  NA{t)  in  number  of  particles  in  wells 
B  and  A  at  any  time  t 

i{NB{t)  -  NA{t))  = 

? “P  (-5&)  (-('VsW  -  ^.4(0)  +  . 

The  polarization  P{t)  due  to  the  applied  electromagnetic  field  is  proportional  to  A^s(t)  - 
A4(t).  By  relating  r  with  2^  exp  and  e*  -  Coo  with  ^  exp  (-^) 

arrives  at  the  Debye  model  (5)  from  atomic  considerations. 
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TliGiG  is  substdntifl,!  icjvsoii  to  believe  tluit  the  behfivior  described  by  the  Debye  model  is 
piessure-dependent.  One  approach  to  understanding  this  pressure-dependence  is  to  extend 
the  above  aiguments  and  consider  the  polarization  from  a  non-eejuilibrium  therinodynamics 
perspective.  A  discussion  of  this  nature  is  given  in  [20]. 

We  however  take  a  different  approach  to  incorporating  pressure-dependence  into  the  Debye 
model.  We  present  the  model  here  and  provide  motivation  in  Section  2.3.  We  begin  by 
assuming  that  the  material-dependent  parameters  in  the  differential  equation  (5)  depend  on 
pressure,  i.e., 


t{p)P  +  P  =  fo(e,(p)  -  eoo{p))E  =  eo(7(p)  -  (:{p))E. 

Wc  suppose  as  a  fiist  aj)proximation  that  each  of  the  pressure-dependent  parameters  can  be 
represented  as  a  mean  value  plus  a  perturbation  that  is  proportional  to  the  pre.ssure 

t{p)  =  To  f  =  To  +  KrP 

lip)  =  7o  +  7  =  7o  +  «7P 

C  (p)  =  Co  +  C  =  Co  + 

Then  the  equation 

t{p)P  -f  P  =  eo(7(p)  -  C(p))^^ 

can  be  written 


(to  -b  Krp)  P  +  P  =  cq  (7o  -  Co  +  (^7  -  K^)p)  E.  (6) 

We  recall  that  the  polarization  term  in  (4)  involves  second-order  time  derivatives.  To  express 
(6)  in  compatible  form,  we  take  the  time  derivative  of  both  sides  to  obtain 

p  -  p  -(  ('yo  -  Co  +  (^7  -  f^<)p)  p  <^o  (s  -  p 

{to  +  K-tP)  (A)  +  Avi-p)  (t'O  +  Av't7?)  ^  ' 

with 


P  =  - 


(ro  -I-  Avp) 


P-f- 


^0  (To  ~  Co  +  (A''-7  -  Ac^)p) 

(To  +  KtP) 


E. 


(8) 


8 


From  here,  we  can  use  (8)  in  (7)  and  then  replace  P  in  (4)  with  the  expression  given  b}'  (7). 

Additionally,  we  should  note  that  the  relation  =  eoo  in  (zi,  1]  becomes  Cr  =  Cooip)  =  Co+f^cP 

in(2i,l]. 


2.3  Motivation  for  pressure-dependence  of  polarization 

The  polarization  described  by  both  the  original  and  pressure-dependent  Debye  models  is  due 
to  the  alignment  of  permanent  dipole  moments  with  the  applied  field.  This  tendency  to  align 
is  inhibited  by  the  presence  of  centrifugal  or  gyroscopic  forces.  These  forces  are  caused  by 
molecular  rotations  and  collisions  [26].  The  pressure  in  a  polarizable  medium  influences  and 
is  influenced  by  the  short  range  particle  interaction  in  the  medium,  including  these  molecular 
rotations  and  collisions.  This  interaction  between  particles  may  serve  to  inhibit  or  facilitate 
the  alignment  of  dipole  moments  with  the  applied  field,  resulting  in  the  modification  of  the 
orientational  polarization  [14].  Figure  3  depicts  this  schematically. 


Figure  3:  Pressure-dependence  of  orientational  polarization 


This  interaction  makes  sense  intuitively;  however,  we  want  to  better  understand  the  interac¬ 
tion  mechanism.  Specifically,  we  want  to  consider  each  polarization  parameter  individually 
and  address  its  possible  pressure-dependence. 
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We  begin  with  the  static  pennitt.ivity  f,.,  which  is  the  electric  permittivity  of  a  dielectric  in 
the  presence  of  a  constant  applied  electric  field.  In  1850  and  1879  res])ectively,  Clausius  and 
Mossitti  independently  discovered  that  for  any  given  material,  the  quantity 


Co 


Mi-') 


is  proportional  to  the  material  density  (page  155,  [18];  page  140,  [25]).  Pressure  variations 
in  the  material  cause  changes  in  its  density.  These  changes  are  reflected  in  the  static  per¬ 
mittivity  of  the  mateiial  due  to  the  law  of  Clausius  and  Mossitti.  So  the  static  permittivity 
can  be  expected  to  depend  on  pressure. 


The  pressure-dependence  of  Cg  does  not  necessarily  suggest  the  pressure-dependence  of  Coo, 
the  permittivity  of  a  material  under  a  very  high  frequency  electric  field.  However  in  the 
inteie.st  of  generalit}',  we  allow  for  the  possibility  that  Is  pressure-dependent.  Pressure- 
independent  behavior  of  Coo  is  just  a  special  case  of  our  model  (see  Section  2.1)  with  =  0 
so  that  Cooip)  =  Co- 

Lastly  we  examine  the  feasibility  of  the  pressure-dependence  of  the  relaxation  parameter 
T.  To  do  so  we  consider  a  dipolar  liquid  which  consists  of  freely  moving  molecules.  If 
an  individual  dipole  changes  its  orientation,  the  nearby  dipoles  shift  to  compensate  and 
produce  a  new  equilibrium  position.  Their  collective  motion  can  be  viewed  as  a  viscous 
frictional  damping  force  that  act.s  on  the  original  dipole.  When  an  electric  field  is  applied, 
its  force  causes  the  dipole  to  align  itself  with  the  field.  The  rate  of  alignment  depends  on  the 
amount  of  friction.  However  since  the  dipole  is  subject  to  the  effects  of  Brownian  motion, 
this  rate  also  depends  on  thermal  fluctuations.  Taking  this  into  account  Debye  derived  the 
following  expression  for  the  relaxation  (page  73,  [2]) 


2kiiT 

where  C  is  the  frictional  constant.  Dipoles  arranged  in  smaller  groups  are  less  apt  to  resist 
reorientation  [13].  This  leads  to  diminished  frictional  effects.  Variations  in  pressure  likely 
alter  the  cohesion  of  dipole  groupings  and  thus  affect  the  friction.  A  specific  example  of 
this  relationship  is  given  in  (page  63,  [21])  for  hard  sphere  fluids.  In  this  case,  the  frictional 
viscosity  constant  is  given  by 


,  _  ksT 
mD 

where  m  is  the  particle  mass  and  D  is  the  self-diffusion  coefficient.  The  self-diffusion  coeffi¬ 
cient  is  pressure-dependent;  that,  is 
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/Trfc^Tx  5 

U  ^  1 

\  m  J 

\pkBT 

J 

where  R  is  the  hard  sphere  diameter,  p  is  the  pressure,  and  p  is  the  liquid  density.  Clearly 
in -this  example  the  relaxation  parameter  r  is  pressure-dependent. 


3  Theoretical  results 


In  order  to  consider  (4)  from  a  theoretical  perspective,  it  is  convenient  to  write  it  in  vari¬ 
ational  form.  We  formulate  the  Gelfand  triple  V  ^  H  ^  V*,  where  H  =  L^(0, 1)  and 
1)  =  {<^  €  H^(0, 1)  :  ^(1)  =  0}.  We  let  <  •,•  >  denote  the  usual  inner 
product.  Then  we  may  write  (4)  as  a  variational  form  of  Maxwell’s  equation  in  second  order 
form 


<  aE{t),  4>  >K.,v  +  <  hE{t),  <j)>  +  <  eP{t),  <f>  > 
+cE{t,0)<i>{0)  +  ay{E{t),cl>)=<F{t),cl>>v,v,  te  [0,r], 
for  all  (j)  G  V .  The  sesquilinear  form  ai  is  defined  by 


(9) 


>, 

where  ^  >  0  is  constant  and  the  parameter  functions  a,b,  and  e  are  determined 

by  the  geometry,  conductivity,  and  instantaneous  polarization  of  the  dielectric.  Since  the 
absorbing  boundary  condition  E  —  cE'  =  0  at  2:  =  0  is  a  natural  condition,  we  incorporate  it 
into  the  variational  formulation  of  the  equation.  The  superconductive  boundary  condition 
at  2:  =  1  however  is  essential  and  is  imposed  in  the  definition  of  V. 

Motivated  by  the  polarization  model  described  in  Section  2,  we  may  rewrite  (9)  as  the  general 
variational  form 


<  aE(t),  (f)  >v',v  T  <  (/)>  +  <  hE(t),  (p  > 

+  <  fo  ■)  ^  > 

+c£(t,0)<^(0)  +  ai(£:(t),<A)  -<  F{t),P>v^y,  0  G  y,  te  [0,T], 

E{0,z)  =  Eoiz)  E(0,z)  =  E^{z), 


(10) 
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where  as  in  (9)  the  coefficients  and  kernel  function  arc  determined  by  system  characteristics. 
For  the  pressure-dependent  Debye  polarization  model  given  in  Section  2,  we  have 

1  +  (too  -  f)/(z,,i)  =  1  +  (Co  + 

1^1  1  ttt(7o-0)  t  (k7  -Ki;)?>((.z))  \ 

I  fo  fn  (.ro+KTp{t,z))  j  ■'(zi.l) 


a{t,z)  = 
hit,  z)  = 


(Here  Iq  is  the  indicator  or  characteristic  function  for  a  set  f^.)  We  note  that  the  sesquilinear 
form  CTi  :  f/  X  K  ->•  R  is  F-continuous  and  F-elliptic,  so  that  there  exist  positive  constants 
Ci,C2  such  that  the  following  inequalities  hold. 


^7l(</>,V0  =  <  9'^',  <  Ci\(/)\y\^h\y 


aiicj),  0)  =  c^  <  0',  0'  >=  c2j0'|2^  >  C2I0I2,. 

In  [9],  the  well-posedness  of  (10)  is  considered  for  general  coefficient,  kernel,  and  forcing 
functions  under  the  following  assumptions 

Al)  The  coefficient  a  along  with  its  derivatives  a  and  d  are  in.  L°°(0,T;  L°°[0, 1]),  and  for  all 
^  €  [0, 1],  a{z)  >  fi(),  for  some  1  >  o-q  >  0. 

A2)  The  coefficient  b  and  its  time  derivative  b  are  in  L°°{0,T:  L°°f0,  ll)  and  b(t  z)  >  0  for 
all(f,z)6[0,r]x[0,l]. 

A3)  The  coefficient  h  is  in  L°“(0,  T;  L“[(),  Ij). 

A4)  The  kernel  function  G  is  in  L“([0,T]  x  [0,T];L“[0, 1]). 

A5)  The  sequilinear  form  is  given  by  cr,  (0, 0)  =  0^  <  0',  0'  >  for  0, 0  G  witfi  c>  0. 
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A6)  The  forcing  function  F  is  in  H^{0,T,V*). 

The  result  is  the  following  theorem.  (See  [9]  for  a  detailed  proof.) 

Theorem  1:  Under  assumptions  Al)-A6),  the  sj''stem  (10)  possesses  a  unique  solution  and 
{E,E)  depends  continuously  on  initial  data  {Eo,Ei)  and  forcing  function  F  from 
iEo,EuF)  eVxHx  H\0,T-,V*)  to  {E,E)  G  L‘^{0,T-,V)  x  L^0,T-H). 

We  may  verify  that  assumptions  Al)-A6)  hold  for  (11)  -  see  [9]. 


3.1  Estimation  of  parameters  in  the  general  variational 
form 

The  previous  well-posedness  result  provides  a  framework  in  which  to  formulate  parameter 
estimation  problems.  As  described  in  Section  1,  the  general  Maxwell  system  treated  by  this 
result  arises  from  a  class  electromagnetic  interrogation  problems.  The  crux  of  these  prob¬ 
lems  is  estimation  of  certain  parameter  values,  namely  dielectric  constants  and  conductivity 
coefficients,  for  the  material  under  interrogation.  The  estimation  problem  typically  involves 
finding  the  parameter  values  that  provide  the  best  fit  between  the  model  and  data  collected 
from  the  actual  system,  using,  for  example,  a  least  squares  criterion.  These  parameter 
estimates  may  then  be  used  to  characterize  the  material. 

In  practice,  the  experimental  data  is  compared  with  finite  dimensional  numerical  approxi¬ 
mations  to  the  model.  In  this  section,  we  examine  the  relationship  between  the  parameter 
estimation  problems  for  the  original  system  (10)  and  for  a  corresponding  finite  dimensional 
system.  We  suppose  that  the  coefficients  and  sesquilinear  form  in  both  (10)  and  its  finite 
dimensional  approximation  depend  on  a  parameter  g  in  a  set  Q.  If  the  exact  solution  to  the 
original  system  (10)  were  accessible,  we  would  consider  the  problem  of  minimizing  the  least 
squares  cost  functional 


Nt 

over  q  ^  Q  where  w  =  is  a  set  of  observations  taken  at  times  tj,  Q  is  &  set  of 

admissible  parameters,  and  O  is  an  observation  operator.  The  form  of  O  depends  on  the 
particular  application  and  set  of  observations.  For  example,  if  Wi  is  a  measurement  of  the 
electric  field  taken  at  a  spatial  point  z  at  time  tj,  then  the  operator  O  entails  evaluations 
of  the  function  E{ti,  qg)  at  a  point  in  space.  Since  we  cannot  obtain  a  closed  form  solution 
to  (10),  we  use  the  solution  E^{t;q)  to  a  finite  dimensional  approximating  system.  The 
solution  E^'{t-,  q)  lies  in  a  finite  dimesional  subset  of  V,  and  satisfies 
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<  a{q)E’^{t),(l)>y.y  +  <  b{(i)E^{t),<t)>  +  <  h{q)E^'{t),cj)> 

+  <!lG{t,sr;q)E^s,-)dsA> 

(13) 

+cE^Mm+^l{(l){E^{t).(l>)  =<  F{t),(l>>y.y 
E^’{Q^z)^V^Eo{z)  E^{0,z)=V^E,{z) 

for  all  <j)  e  V^.  In  particular,  we  define  the  piecewise  linear  basis  elements  with 

nodal  values  (l>f{k/N)  =  4;,  fc  =  0, 1, . . . ,  Af,  and  let  =  span  c  K 

Then  we  define  to  be  the  quasi-L^O,  1)  projection  (see  [23],  [4],  [12])  of  V*  onto  V^' 
defined  by 


<  V^v\v^  >N=<  V*,  >y.y  for  V*  G  K*  and  for  all  e 


where 


<»",«">«=  /' /"(«,  V)(z) 

Jo 

and  is  the  nodal  value  linear  interpolation  operator  for  .  It  is  shown  in  [23]  that  the 
operator  "P^  is  well-defined  and  satisfies 


<  K'Mh  for  <I>  ell 

\V^(l)\y  <  iq\(j)\v  ior  (/)  e  V. 


As  expected,  the  corresponding  cost  functional  for  the  finite  dimensional  system  is 


Nl 

J^(r/,  7/;)  =  \OE^{ti]  q)  -  Wif  .  (15) 

i-l 

Again  the  form  of  the  operator  O  is  chosen  to  correspond  to  the  type  of  data  collected. 

In  [9],  we  established  the  well-po.sedness  of  (10)  with  .solutions  E  in  L-{0,T;V)  and  E  E 
F,  H),  where  Y  —  //)j(0, 1)  and  H  =  /y^(0, 1),  and  we  also  verified  that,  a  unique 
solution  to  (13)  exists.  These  results  hold  provided  that,  for  each  q  E  Q,  Assumptions  Al)- 
A6)  are  satisfied  and  the  sesquiliuear  form  a\  is  V-continuous  and  V"-elliptic.  Moreover  in 
[24],  we  .show  that  the  solution  E  of  (10)  lias  the  enhanced  regularity  E  E  H'\0,T;V*)  n 
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//■^(0,T;if)  n  H\0,T-,V)  under  consistency  conditions  for  the  initial  conditions  and  the 
following  assumptions: 

A7)  The  second  time  derivative,  b,  of  6is  in  L°°(0,T;T~[0,1]). 

A8)  The  first  and  second  time  derivatives,  h  and  h,  of  h  are  in  jL°°(0,  T;  L°°[0,  1]). 

A9)  The  first  and  second  derivatives  with  respect  to  the  first  temporal  variable,  and 

^G,  of  the  kernel  function  G  are  in  L°°((0,T)  x  (0,r);  L°°[0, 1]). 

AlO)  The  forcing  function  F  is  in  H'^{0,T,V*)  and  is  of  the  form  F{t,z)  =  git)5{z)  with 
g{t)  e  and  ^(0)  =  ^(0)  =  0.  (This  assumption  replaces  A6).) 

(Verification  of  these  assumptions  A7)-A9)  for  the  pressure-dependent  Debye  polarization 
model  are  also  given  in  [24].) 

We  now  make  the  following  assumptions  about  the  set  of  admissible  parameters  Q,  the  state 
space  and  the  projection  operator  . 

Bl)  The  finite  dimensional  set  Q  lies  in  a  metric  space  Q  with  a  metric  d  and  is  compact 
with  respect  to  this  metric. 

B2)  The  finite  dimensional  subspaces  are  subsets  of  V. 

B3)  For  each  (j)  G  V,  \(j)  —  — y  Q  ds  N  — y  oo. 

B4)  For  each  (/>£  (})\h  ->  0  as  ->  oo. 

Verifications  of  B3)  and  B4)  for  our  particular  are  given  in  [23],  We  now  make  a  further 
assumption  on  the  sesquilinear  form  ui.  We  assume  that  ai  =  ai(g)  is  defined  on  Q  and 
satisfies 

HI) 

ki(9i)(<A.'0)  -cfi(<j2)(A </')!  <  'yi(<iu‘i7)\4>\v\’l>\v 
for  qi,q2  E  Q  where  7  depends  only  on  Q. 

For  the  electromagnetic  system  in  consideration  in  this  paper,  Assumption  HI)  is  unnecessary 
since  ai  is  independent  of  q.  However  for  the  purpose  of  establishing  a  more  general  result, 
we  do  not  assume  here  that  our  sesquilinear  form  is  parameter  independent. 

Furthermore  we  make  the  following  assumption  about  our  coefficients. 

All)  The  coefficients  depend  continuously  on  q  so  that  as  d{q,q^)  -4-  0,  we  have 
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i)  |o(q)  -  fl((7")|/,<x. ->  0 

n)  \h{q)-b{q^)\r^^^() 

Hi)  \h{q)  -  h{q^)\,^o. 

w)  \G{q)-G{q^)\,,^-^Q. 

The  above  continuity  along  with  the  compactness  of  Q  implies  that  the  images  a{Q),  b{Q), 
h{Q),  and  G{Q)  are  compact.  Thus  each  coefficient  can  be  bounded  independently  of  q.  We 
assume  throughout  that  all  bounds  on  our  coefficients  do  not  depend  on  q. 

By  solving  the  parameter  estimation  problems  related  to  (13),  (15)  we  obtain  a  sequence  of 
estimates  {g  }.  We  wish  to  demonstrate  that  under  certain  conditions  this  sequence  (or  a 
subsequence)  converges  to  the  estimate  corresponding  to  the  problem  related  to  (10),  (12). 
In  order  to  do  this,  we  .state  the  following  claim,  which  can  be  found  (along  with  a  proof) 
as  Theorem  5.1  in  [11]. 

Theorem  2:  To  obtain  convergence  of  at  least  a  subsequence  of  to  a  solution  q  of 
minimizing  (12)  subject  to  (10),  it  .suffices,  under  assumption  Bl),  to  argue  that  for  arbitrar}^ 
sequences  {q^}  in  Q  with  q^  ->  q  in  Q,  we  have 


OE^if,q^)-^OE{t;q). 


In  [11],  the  operator  O  is  general  enough  to  include  functions  that  map  functions  /  such 
that  /  :  T  -4  V  to  the  space  of  observations  W,  where  T  is  an  appropriately  chosen  (see  [11] 
and  [8])  subset  of  [0,T]  that  contains  the  times  of  observation  and  V  is  a  space  containing 
E{t,  •).  In  the  numerical  examples  presented  in  this  paper,  the  observations  correspond  to 
the  values  of  the  electric  field  at  the  point  ^  =  0  at  various  times,  i.e.,  {£'(ij,0)};  thus  the 

operator  O  involves  pointwise  evaluation  of  E  at  many  points  in  time  and  one  specific  point 
in  space. 

We  suppose  that  and  satisfy  B2)-B4),  the  sesquilinear  form  a,  satisfies  Hi),  the 
coefficients  satisfy  assumptions  Al)-All)  and  we  let  q^  €  Q  be  arbitary  such  that  q^  — >  q 
in  Q.  Our  primary  goal  is  to  show  that  as  A'  -)■  00 


E^'{t,0-,g^')-^E{t,0-q)  (16) 

for  each  t  e  [0,T].  However,  hero  we  verify  a  more  general  result.  We  show  that  for  each 

f,G[o,r]  ' 
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(17) 


E^(t]q^)  — >  E{t-,q)  in  the  V  norm 
E^{t-,  q^)  — >  E{t]  q)  in  the  H  norm 

as  N  oo,  where  E^,E^  are  the  solutions  to  (13)  and  E,E  are  the  solutions  to  (10).  We 
note  that  we  may  evaluate  these  functions  pointwise  in  t  due  to  the  enhanced  regularity  of 
solutions.  Moveover,  using  the  equivalence  of  norms,  we  see  that  (17)  implies  (16)  and  we 
have  the  result  we  need  for  our  computations.  . 

We  point  out  that  for  a  sequence  q^  ~  q  for  all  N,  the  desired  result  implies  convergence  of 
the  finite  dimensional  approximation  to  the  true  solution.  This  is  important  when  considering 
numerical  approximations  to  the  solution. 

We  have  established  previously  that  the  solution  of  (10)  satisfies  E{t)  e  V  and  E{t)  e  H  for 
each  t.  Since 

.\E^{t-,q^)-E{P,q)\v  <  \E^ E{t-,q)\v  +  \V^ E{t;q)  -  E{t-,q)\y 
and  B3)  guarantees  \E^E{t',  q)  —  E{t\  q)\v  — >  0  as  oo,  we  need  onl}'’  show  that 

\E^it;  q^)  -  V^E{t;  q)\v 0  as  N  oo 
for  each  t  E  [0,  T].  In  the  same  way,  it  suffices  to  show  that 

\E^it;q^)-V^E{t-q)\H-^OasN-^oo 
for  each  t  E  [0,  T]  to  obtain  the  second  result. 

We  let  E^  =  E^{t;  q^),E  =  E{t;  q),  and  =  E^{t]  q^)  -  E{t]  q). 

Subtracting  (10)  from  (13),  we  have  for  4>  eV^ 

<  a{q^)E^  -  a{q)E,  <i>  >  +  <  b{q^)E^  -  b{q)E,  (j>  > 

+  <h{q^)E^  -h{q)E,(l)> 

+  <  fo(G(q^)E^  -  G(q)E)  ds,  >  +c{E^{t,  0)  -  E{t,  0))<^(0) 
+ai(g^)(S^^)-ai(g)(E,</.)  =  0. 
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We  add  and  subtract  E  and  its  derivatives  and  rearrange  terms  to  obtain 


<  a{q^’){E^  -V^E),cj,>  -  V^E,  cf>) 

-<  a{q)E,^>  -  <  a{q’^)V^E^(f)>  +a^{q){E,(j))  -  ai{q^){V’^ E,  (/)) 
+c{Eit,  0)  -  V^E{t,  0))(j>{Q)+  <  Kq)E  -  b{q^)E^,  (j>  > 

+  <  h{q)E  —  h{q^)E^ ,(j)  >  +  <  G{q)E  —  G{q^)E^'  ds^cj)  >  . 

We  choose  the  test  function  </>  =  A^'  G  so  that 

<  a{q^){E^  -V^ E),K^  >  +o,{q^^){E^  -V^E,K^) 

+c{E^{t,  0)  -  T’^E{t,  0))A^(f,  0)  = 

<  a('7)^,A^  >  -  <  a{q^')V^E,k^  >  +Oi{q){E,K^)  -  ax{q^){V^ E,  K^) 

+c{E{U  0)  -  V^E{t,  0))  A^(/,,  0)+  <  b{q)E  ~  b{q^)E^,  > 

+  <  h{q)E  -  h{q^)E^,  >  +  <  fl,iG{q)E  -  G{q^')E'^)  ds,  A^'  >  . 

We  note  that 

2  <  a{q‘^'){E^  -  V'^E),A^  >  +2ariq^){E^  -  V'^'E,A^') 

=  i{<  o(9'^)A'",  A"  >  +e7i(e/^)(A'^,A^))  -  <  o(e7^^)A^,A^  >  . 

Tlien  we  have 
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ki  (<  >  +a:(?^)(A^  A^))  +c(A^(f,0))2  = 

<  a{q)E,  A^  >  -  <  a{q^)V^E,  >  +  <  a{q^)K^ ,  > 

+ar{q){E,L^)  -  a,{q^){V^ E, 

+c(E(f,  0)  -  V^E{t,  0)) A^(^,  0)+  <  b{q)E  -  b{q^)E^,  A^  > 

+  <  h{q)E  -  h{q^)E^,  A^  >  +  <  fo(G(q)E  -  G(q^)E^)  ds,  A^  >  . 
Integration  with  respect  to  t  yields 

<  a(g^)A^(t),A^(t)  >  +ai(g^)(A^(t),  A^(^)) 

+2/;c(A^(e,0))2rfe  = 

2  /o  {  <  Ci{q)E,  A^  >  -  <  a{q^)V^E,  A^'  >  +  <  a(g^)A^,  A^  > 

+ai{q){E,A^) -ai{q^){P^E,A^) 

+c(S(e,  0)  -  V^E{^,  0))A^(e,  0)+  <  b{q)E  -  b{q^)E^,  A^  > 

+  <  h{q)E  -  h{q^)E^,  A^  >  +  <  //  G{q)E  -  G{q^)E^  ds,  A^  >  } 

+  <  a(5^)A^(0),A^(0)  >  +c7i(9^)(A^(0),  A^(0)). 

We  now  use  the  definition  of  A^  to  obtain 

A^(0)  =  £;^(0)  -  P^E{0)  =  ^^(0)  -  V^Eo  =  0 

A^(0)  =  A"(0)  -  V^E{0)  =  E^{0)  -  V^Ei  =  0. 

We  may  then  write 
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+2/„‘c(A^(e,0))2f/e  = 

2  {  <  a{q)E,  A"  >  -  <  a{q^)V^E,  A^  >  +  <  a{q^)A^,  A^  > 

+<Ti  {q){E,  A^)  -  ai  {q^){V^E,  A^) 

+c(E(^,  0)  -  V^Ei^,  0))A^(^,  0)+  <  b{q)E  -  b{q^)E^,  A"  > 

+  <  h{q)E  -  h{q^)E^,  A^  >  +  <  G{q)E  -  G{q^)E^  ds,  A^  >  }  d^. 

In  order  to  bound  the  right  side  of  (18),  we  derive  tlie  following  estimates: 

Estimate  1: 

/o  (2  <  a{q)E,A^  >  -2  <  aiq^)V^E,A^  >  +  <  a{q^)A^,A^  >  )  d^ 

=  fo  (2  <  («(<?)  -  a(q^))B,  A^'  >  +2  <  a((/^)(E  -  P^E),  A"  > 

+  <  d{q’^)A^,  A"  >  j  d^ 

<  fo  l(«(9)  -  +  1(5  +  |a(9^)|2,)|A^'i|,  d^ 

<  Kq)  -  «('^^)lioo  j„'  jEll  d^  +  |a((/^)|2^  jE  -  V’^E\l  d^ 

+  |(5  +  |%^)lioo)/„^|A^||,r/e. 

Estimate  2: 

2f^^ay{q){E,A^)-a,{q^[){V^’E,A^)  d^ 

=  2ioai(g")(P"E,A^') -f^i(<7)(^,A^)  d^ +  2{a,{q){E{t),A^ {t)) 

-a,{q^’){V^E{t\A^{t))) 
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=  2S‘ar{q''){P''E  -  E,A«)+^,(q«)(B,A^)  -  <7,(?)(B,  A")  d( 

+2^<7i(q)(E(t),  A''(t))  -  a,  (q'')(E(t),  A''(t)) 

+<7,(?'')(£(t) -■?"£(().  A''(())) 

<  cfiP-'i;  -Bil-h  7Hd(q,  q^nWv  +  aiA"!?,  d( 

+^lP''B(t)  -  E(t)ll  +  f  (d'h,  ?''))^|B(i)Pv  +  2e|A''(f)Pv. 

where  e  >  0  is  arbitrary. 

Estimate  3: 


2cm,0)-r^B(^,0))A^(^,0) 

<  c2|E(e,  0)  -  0)P  +  I  A^(e,  0)  p 

(Here  we  use  the  fact  that  |^|f/  is  equivalent  to  |</)'||^  +  |^(0)p  so  that  \(l)\y  >  K{\(j)'W  + 

i.ji(o)n>/?i^(o)p.) 

Estimate  4: 

2j^<b{q)E-b{q^)E^,A^> 

=  2  /o‘{<  b{q)iE  -  r^E),  >  +  <  ib{q)  -  b{q^))V^E,  A^  > 

+  <  6(9^)A^,A^  >} 

<  |6(9)(E  -  V^E)\\  +  +  |(6(g)  -  b{q^))V^E\\  +  3|A^|2,  d^ 

<  m\U  fo  1^  -  +  3)  fo' \A^\l 
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Estimate  5: 


2/;  <%)£-%")/;»,  A"  >  rff 
=  2/„‘{<  /i(9)(£; A"  >  +  <  A'''  > 

+  </ib")A'',A''  >) 

<  lo  (l/‘M(E  -  +  |/i(,,»)A«||,  +  slA"!!, 

+I(A(';)  -  A(«''))'P''B||,)  (if 

<  IhWIi- 1^; -  ^"4.11,  <'5  + IM7'')II-/.‘|a''|5,  de  +  3/‘|A''||,  d? 

+|A{5)  -  /.fa'')|i„  /„■  df. 

Estimate  6: 

2  Jo  <  //(G(d)B  -  G(,'*)B'')  d.5,  A"  > 

=  2  /;  <  /«  GW(B  -  P^B)  ds,  A"  >  +  <  C„«(GW  -  G(,''))P'«£  do,  A"  > 

+  <//G(?'')A'''ds,A''>  df 

<  /„'{l  Jo*  G(,)(iS  -  P^B)  d,g  +  I  /„«(G(,)  -  G(,/'»))P«f;  d,,[f, 

+|/„«G(,'')A''d,,|f,  +  3|A''|i,)d{ 

=  /o'  /o'  I  f«  G(?)(B  -  P"B)  doP  do  d? 

+  /) /o'  I  loi^io)  -  G{q^))V^E  o/sp  dz 

+  /o /o'  I  /<fG(d'')A«  d,,p  dz  df  +  3/;  lA"!?,  oif 

<  |G{?)|i«  I  j„*(B  -  p"/;)  dsp  dz  d(  +  |G(<7)  -  G(<7'')P^  fo  Id  I  Io 'P" E  *1’  *  li? 

+|G(d'')li./,7o'  l/jA"  dz  df +  3./;;|A"||,  d^ 
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<  \G(q)\l„Tf‘  |B  -  dz  di  +  |G(5)  -  G(?'')||^7’/„‘ /J  dz  d( 

<  lG{q)iUT^f‘lB-V''E\ld(  +  \G{q)-G(q'<)\l„Pf‘\P^E\%d^ 
+|G(9'')I!»tV„‘|a''|5,  de  +  3/<;|A'»i5,  di. 

Using  these  estimates,  Assumption  HI),  the  U-continuity  and  V -ellipticity  of  cti,  and  the 
fact  that  <  \(j)\y,  we  may  rewrite  (18)  as 

|^a(g^)A^(t)|2^  +  C2|A^(0|f/  +  2  /oC|A^(e,  0)p 
<  5f(^)+5|'(i)  +  2e|A^(i)|?. 

+  /o  {lA^fco)!"  +  (f  +  +  l|d(,»)|i„)|A''|l, 

+  (2  +  l'‘(9'')li«  +  nG(q^)\l«)  |A»P,,}  di 

where 

df'W  =  f‘{\a(q”)\UE-T^E\l 

+(c;  +  cW,)|B  -  +  |K9)|i_|B  -  P’^Ell 

+  (|M9)li»  +  J’"|G(9)li»)  \E  -  df 


=  Hq)  -  a{g'')\l^  \E\l  ,li 
Jo 

+7"fl(<7,9^)  [  \E\l  +  |6(fy)  -  b{q^)\l^  I  \V^E\l 
Jo  ./o 

+  Ol‘M  -  '•(«'')li»  +  T^IG(q)  -  G(q’')ll^]  x 

/'lP’'Kdi+^d(q,q'')lBll)ll 

Jo  ^ 

Since  e  >  0  is  arbitrary,  we  may  choose  it  to  be  sucli  that  1  >  C2  -  26  >  0.  Furthermore,  the 
wave  speed  c  satisfies  2c  >>  1.  We  then  use  Assumptions  Al)-A4)  to  claim  that  there  exist 
constants  />'i,i/2  >  1  and  1  >  ao  >  0,  independent  of  q,  such  that 


«o|A^(t)|2,  +  (C2  -26)|A^(t)|f,  + /^*(2c-  l)iA^(e,0)|2  d^ 

Finally  recalling  the  bounds  on  1^2,  (^o,  and  C2  —  26,  we  may  rewrite  the  ineciuality  as 

ao|A^(t)|2,  +  (c2-26)|A^(/,)|2 

<  «  +  j'  (^)  r, lA"!?,  +  d( 

oolA^I;,  +  (c2  —  2€)|A''  \y  df;. 


<  S^'il)  +  S^(l.)  + 


(h{('2  ~  26) 


In  order  to  apply  Gronwall’s  inequality  to  obtain  uniform  convergence  of  A^  and  A^  in  t  as 
N  00,  we  must  establish  the  uniform  convergence  of  and  .  We  have  from  B3)  and 
B4)  that  \E{t)  —  E{t)\ji,  \E(t)  —  V^E(t)\Y,  and  |£'(t)  —  V^Elf^ly  converge  to  zero  as 
Af  -»  c»  for  each  /.  Since  this  convergence  is  dominated  and  {A'(t)}te[o,7’]  is  compact  in  V, 
we  have  that 


->  0  uniformly  in  t  rs  N  00. 

Moieovei,  the  boundedness  of  E,  E,  anrl  E  given  by  the  enhanced  regularity  results  and  the 
assumption  AS)  im])ly  that 
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->  0  uniformly  in  i  as  ->  oo  and  q  in  Q. 

Then  vve  may  apply  Gronwall’s  inequality  to  conclude  that 

^^Pfe[o,T]  l^^(^')lv  0  as  A/^  — >  oo 

suPt£[o,T]l^^(Olff  OasN-^oo 

which  is  sufficient  to  prove  the  desired  result. 


3.2  Estimation  of  parameters  in  the  system  with  pressure-dependent 
Debye  polarization 

The  geneial  system  (10)  is  formulated  to  accomodate  systems  arising  from  a  variety  of 
electromagnetic  interrogation  problems.  We  are  concerned  here  with  a  particular  system 
that  incorporates  a  pressure-dependent  model  for  Debye  polarization.  We  demonstrate  that 
this  system  satisfies  Assumptions  All),  Bl)-B3)  and  HI)  and  thus  that  the  results  of  the 

previous  section  apply.  (We  note  that  verifications  of  Assumptions  A1)-A10)  are  given  in  [9l 
and  [24].)  /  ^  6  l  J 

The  system  we  wish  to  consider  is  given  by  (10)  wdth  the  parameter-dependent  coefficients, 
kernel  and  forcing  functions,  and  sesquilinear  form  (11). 

For  this  system,  the  set  of  admissible  parameters  Q  a  subset  of  R^,  where  seven  is  the 
number  of  parameters  to  be  estimated  (in  addition  to  the  six  polarization  parameters  from 
Section  2,  one  is  often  interested  in  estimating  the  conductivity  coefficient  a).  Here  we 
consider  g  G  Q  C  where  q  =  (o',  70,  Tq, /v.^, /c^,  k^).  We  choose  the  admissible  set  Q  to 
insuie  that  our  Debye  coefficients  are  well-defined.  First,  because  of  the  physical  meaning  of 
these  parameters,  the  values  of  a,  70,  Co  and  tq  must  be  positive.  Then  for  a  given  pressure 
\vave  p  with  p,p  E  L°°{0,  T;  T°°(0, 1))  and  a  fixed  5  >  0,  we  admit  only  values  of  Kj,  k^,  and 
K-r  such  that  7o  +  K^p{t,  z),  Co  +  z),  and  tq  +  Hrp{t,  z)  are  greater  than  S  for  all  2  G  [0, 1] 

and  t  G  [0,r].  In  addition  to  these  requirements,  we  assume  that  the  admissible  parameter 
set  is  closed  and  bounded  in  R'^. 

We  recall  that  q  ^  in  the  standard  Euclidean  metric  is  equivalent  to  the  convergence 
of  each  component  of  q^ .  Moreover,  any  closed  and  bounded  sets  Q  in  R^  are  compact 
and  satisfy  the  conditions  of  Bl).  The  conditions  B2)-B4)  are  satisfied  by  ,  which  is 
in  this  case  the  set  of  finite  dimensional  linear  piecewise  basis  elements,  and  the  projection 
operator  .  To  verify  Hi),  we  note  that  cFi{q){(j>-,i))  =  <  (p'yijj'  >  is  independent  of  q 

and  \(7i(qi){(l),ip)  -  ai{q2){(j),'ip)\  =  0  for  any  91,92  G  Q. 
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We  next  verif}'  All)i)-iv)  for  the  coefficients  in  our  model.  We  note  that  as  — >  q  we  have, 

for  a  given  e.  L°°, 


\a^  -  (7|  ->  0 


(19) 


Co  P  fCo  +  0 


(20) 


7o^ 


Co  +  f'x)p  “  ^7o  ~  Co  +  {fty  — 


0 


(21) 


+  >^tP  -  (jO  ~ 


^0 


L°= 


(22) 


(k"  ~  K^)p  -  ((k,  -  K^)p'j 


loo 


(23) 


and 


(24) 


We  use  (20)  directly  to  claim  that  In)?")  -  a((;)|,,«  ->  0  whenever  q'^  -t  q;  thus  AU)i) 
holds. 

In  demonstrating  that  All)ii)  holds,  we  observe  that 


IM?")  -  K<7)I 


< 


^0 


.A' 


a 


+  sup  sup 

26(21,1] 


7^  Co  +  (^7  -  ^f)p(f  ,^)\  _  /  7o  -  Co  +  (Ky  -  K^)p{t,z) 

n-  A^  rr\(  -t-  I  1  _  1  /  I  \ 


+  '^^TP(f>  '2) 


To  +  KrP(t,  Z) 


Then  we  may  apply  equations  (19),  (21),  and  (22)  and  the  quotient  rule  of  limits  to  conclude 
that  |/>((i^)  —  b{q)\ioo  — ^  0  as  q^  — >  q  and  All)ii)  is  satisfied. 

BVom  the  definition  of  h,  we  have 


26 


|/i(g^)  —  h(q)\L^  <  sup  sup 

«e[0,T]  Z6[zi,l] 


-  K^)p{t,  z) 
+  K^p{t,  z) 


(Kj  -  K^)pit,  z) 
To  +  Hrp{t,  z) 


(1  +  Krpjt,  Z))(0'0  -  Co  +  (Kj  -  K^)p{t,  z)) 

(to  +  KrP{t,  z))2 


-  C  + 

K^p{t,z)) 


From  equations  (23),  (22),  (24),  and  (21)  and  the  product  and  quotient  rules  of  limits,  we 
may  conclude  that  All)iii)  holds. 

To  show  that  All)iv)  holds,  we  argue  that 


~  G{q)\Lcc  < 


sup  sup  exp 
(«,«)e[o,r]x[o,r]ze[zi,i] 


-d^  \i'^  +  i^rPis,  z)){jo  -  C  +  )p(s,  z)) 

To^  +  <p{^,  z)  )  +  <p(i,  z)Y{t^  +  Kfp(s,  ^))2 


\  (1  +  KrPjs,  ^))(7o  -  Co  +  (/C7  -  /t<)p(s,  z)) 

,  To  K.rP{(,,z)  j  {Tq  + K,TP{t,z)Y{To-\- K.rP{s,z)Y 


We  note  that  equation  (22)  coupled  with  the  quotient  rule  for  limits  allows  us  to  assert  that 


T^  +  k!^p{^,z) 


-dj 

To  +  KtP(C,  z) 


^  Q-  Thus,  we  may  use  equations  (24),  (21),  (22),  and  (25)  with  the  quotient  and 

product  rules  for  limits  to  verify  that 


|G(g^)-G(g)U«->0 

as  q^  q. 

We  have  thus  verified  that  the  theory  established  in  Section  3.1  can  be  applied  to  the  system 
corresponding  to  pressure-dependent  Debj'^e  polarization  described  by  the  coefficients  (11). 
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4  Numerical  examples 


In  the  previous  section,  we  demonstrated  that  a  (sub)sequence  of  minimizers  {q^}  of  the  cost 
functionals  (15)  converges  to  a  mininiizer  g  of  (12).  In  this  section  we  present  computational 
results  for  the  problem  of  finding  ry"  for  a  fixed  N.  (We  attempt  to  choose  AT  large  so  that 
is  close  to  q.) 

This  problem  is  equivalent  to  our  main  objective,  estimating  the  polarization  and  conductiv¬ 
ity  parameters  of  a  dielectric  by  comparing  numerical  solutions  of  the  model  with  experimen¬ 
tal  data.  We  recall  from  Section  2  that  our  polarization  model  has  six  material-dependent 
parameters;  in  this  section  we  fix  a  and  consider  only  these  six  variables  tq,  Co,  7o,  and 

in  the  equation 


(to  +  K-^p)  P  -f  P  —  fo  (70  —  Co  +  («^7  —  K^)p)  E. 

We  want  to  test  the  feasibilty  of  estimating  them  from  experimental  data.  At  this  time,  we  do 
not  yet  have  data  from  experiments  (an  experimental  device  to  obtain  such  data  is  currently 
being  constructed);  instead  we  create  simulated  data  from  our  computations.  The  simulated 
data  consists  of  the  boundary  data  from  a  numerical  airproximation  to  the  solution  of  the 
S3'Stem  with  added  noise.  (See  [9]  for  sample  numerical  solutions  to  the  forward  problem.) 
We  compute  this  approximate  solution  with  fixed  parameters  values.  These  values  are  then 
thought  of  as  our  “unknown”  true  material  parameters.  The  goal  is  to  estimate  these  values. 
We  appraise  our  ability  to  solve  the  problem  by  comparing  the  e.stimates  wdth  the  true  fixed 
values.  If  we  cannot  accurately  approximate  the  parameter  values  in  this  context,  we  cannot 
expect  to  be  able  to  estimate  them  in  an  experimental  setting. 

We  let  q  generically  denote  the  set  of  i)aranicters  we  wish  to  estimate  in  the  examples 
presented  below;  these  may  include  the  mean  values  in  the  polarization  model,  tq,  70,  and 
Co  and/or  the  coefficients  of  pressure  in  the  polarization  model,  Kr,  and  We  let  q* 
denote  the  true  values  of  the  corrsponding  “unknown”  parameters.  We  leave  the  values  of 
all  other  parameters  fixed. 

There  are  two  sets  of  electromagnetic  reflections  that  reach  the  boundary.  The  first,  after 
the  initial  signal,  are  tlie  reflections  from  the  air/dielectric  interface  and  the  second  are  from 
the  viitual  interface  produced  by  the  acoustic  pressure  wave.  (Figure  4  depicts  each  set  of 
reflections  separately.)  In  some  scenarios,  using  data  that  contains  only  one  set  of  reflections 
may  be  advantageous.  For  example,  one  may  use  the  data  from  the  initial  signal  and  the 
reflections  from  the  air/dielectric  interface  (i.e.,  the  first  section  of  data  in  Figure  4)  to 
refine  the  initial  parameter  estimates  and  then  use  these  refinements  with  the  data  from  the 
acoustic  interface  reflections  to  obtain  final  estimates.  In  another  approach,  one  may  use 
just  the  data  from  the  acoustic  interface  reflections  l,o  estimate  the  parameters.  In  any  case, 
the  cost  functional  for  t,he  examjrles  given  here  is  of  the  form 
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iei 

where  /  corresponds  to  an  appropriately  chosen  data  set.  (Since  here  we  consider  exclusively 
the  finite  dimensional  system  for  a  fixed  N,  we  drop  the  N  for  ease  of  notation.) 


Figure  4:  The  two  sets  of  reflections  that  reach  the  boundary 


We  used  a  Nelder  Mead  optimization  routine  [19]  to  find  the  parameter  values  that  minimize 
the  cost  function.  This  optimization  method  is  a  gradient-free,  simplex  search  method.  The 
use  of  a  gradient-based  method  to  solve  this  problem  is  impractical  due  to  computational 
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time.  The  particular  routine  used  requires  an  initial  simplex  of  estimates  and  a  termination 
tolerance  for  the  difference  between  subsequent  function  evaluations.  We  choose  the  initial 
estimates  to  have  varying  levels  of  error  in  relation  to  the  true  parameter  values.  For  these 
computations,  we  set  tlic  termination  tolerance  at  le-09. 

As  already  noted,  we  created  simulated  data  to  test  our  algorithms.  The  data  set  without 
noise  is  simply  observations  at  the  boundary  taken  from  a  forward  simulation  of  the  model 
using  the  parameter  set  q* .  The  data  sets  with  error  were  created  by  adding  an  appropriate 
amount  of  normally  distributed  relative  random  noise  to  the  original  data  set.  The  random 
noise  was  generated  by  the  MATLAB  command  randn  which  creates  normallj'^  distributed 
noise  with  mean  0  and  variance  1  and  was  scaled  and  shifted  appropriatel}^  Because  the 
noise  is  relative,  the  magnitude  of  noise  is  greater  in  the  intervals  of  data  that  contain  the 
initial  interrogating  impulse  and  the  reflections. 

We  next  present  sample  results  for  specific  parameter  estimation  problems.  We  first  consider 
the  problem  of  estimating  q*  =  [70*, Co,  A*  =  -  [78.2,5.5,0.10545728042059]  from 

data  with  varying  levels  of  noise.  Here  Tq  is  so  small  that  it  is  advantageous  to  estimate  A*, 
a  scaled  function  of  Tq  ;  an  estimated  value  of  Tq  may  be  computed  from  an  estimation  of 
A*.  We  use  the  data  containing  the  initial  signal  and  the  reflections  from  the  air/dielectric 
interface  to  refine  the  initial  parameter  estimates  <70  and  the  data  containing  the  reflections 
from  the  acoustic  interface  to  obtain  final  estimates.  We  present  the  results  in  the  following 
table. 
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Initial  estimate  go 

Final  estimate  for  data  without  noise  g 

go  =  l.Og* 

[78.2,  5.5,  0.10545728042059] 

go  =  0.95g* 

[78.190631,  5.499999,  0.105462] 

go  =  1.05g* 

[78.197841,  5.500000,  0.105458] 

go  =  0.9g* 

[77.128609,  5.499937,  0.107040] 

go  =  l.lg* 

[70.763992,  5.499504,  0.117540] 

Final  estimate  for  data  with  1%  noise 

go  =  0.95g* 

[78.210472,  5.499997,  0.105444] 

go  =  1.05g* 

[77.260485,  5.499940,  0.106852] 

go  =  0.9g* 

[78.198944,  5.490998,  0.105461] 

go  =  l.lg* 

[78.482467,  5.500010,  0.105052] 

Final  estimate  for  data  with  5%  noise 

go  =  0.95g* 

[74.876211,  5.499764,  0.110559] 

go  =  1.05g* 

[77.975797,  5.499939,  0.105722] 

go  =  0.9g* 

[78.337467,  5.500010,  0.105186] 

go  =  l.lg* 

[78.405972,  5.499987,  0.105175] 

Table  1:  Parameter  estimation  results  for 

9*  =  =  [78.2,5.5,0.10545728042059] 


These  results  illustrate  that  it  is  possible  to  recover  accurate  approximations  of  7q,  ^q,  and  A* 
in  the  presence  of  noise  and  with  error  up  to  10%  in  the  initial  estimates.  A  few  of  the  results 
are  unexpected,  for  instance  the  ability  to  approximate  the  values  better  in  the  presence  of 
5%  noise  with  an  initial  guess  with  -10%  error  than  with  an  initial  guess  with  -5%  error.  We 
suspect  these  anomalies  are  due  to  the  simplex  search  nature  of  the  optimization  routine. 

The  results  in  the  tables  clearly  indicate  that  we  can  recover  q*  without  much  error.  However 
it  is  often  illustrative  to  compare  the  solutions  calculated  with  the  estimates  with  the  solu¬ 
tions  calculated  with  q*.  To  do  this,  we  plot  the  absolute  value  of  the  error  for  the  solutions 
computed  at  the  boundary,  i.e.. 


where  q  is  the  final  estimate  (given  in  the  table).  As  an  example.  Figure  5  depicts  this 
error  for  the  estimation  problem  with  5%  noise  and  an  initial  guess  with  -10%  error.  We  see 
that  overall  the  magnitude  of  error  is  small  and  that,  as  expected,  the  most  error  occurs  in 
approximating  the  material  and  acoustic  interface  reflections. 
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Figure  5:  |F(fi,0;9)  -  E{ti,Q-,(j*)\  vs  U  - 
Absolute  error  for  the  parameter  estimation  problem 
with  5%  noise  and  an  initial  guess  with  -10%  error. 


VVe  next  consider  the  estimation  of  q*  =  [«;,  «*,  ac;]  =  [46.92, 1.65, 1.581 139e  -  09]  from  the 
previous  data  sets.  Since  these  parameters  are  the  coefficients  of  pressure,  they  are  irrelevant 
and  undeterminable  until  the  electromagnetic/acoustic  interaction  occurs.  Thus  we  include 
onljf  the  second  section  of  data  in  the  cost  functional.  We  present  the  results  below. 


Initial  estimate  q^ 

Final  estimate  for  data  without  noise  q 

qo  =  0.99g* 

[46.450800,  1.633500,  3.16174e-09] 

qo  =  l-Olq* 

[46.92,  1.65,  1.58114e-09] 

qo  =  0.95q* 

[  44.574000,  1.567500,  5.52700e-09] 

qo  =  1.05g*  . 

[49.266000,  1.732500,  0] 

qo  =  0.9q* 

[42.228000,  1.485000,  0] 

qo  =  1-1(7* 

[51.61200,  1.815000,  0] 

Final  estimate  for  data  with  1%  noise 

qo  =  0.99g* 

[46.450800,  1.633500,  3.17078e-09] 

qo  =  1-Olg* 

[47.389200,  1.666500,  -7.2435e-10] 

qo  =  0.95(7* 

[44.574000,  1.567500,  5.53493e-09] 

qo  -  l-05g* 

[49.266000,  1.732500,  -1.064074e-08] 

qo  =  0.9q* 

[42.228000,  1.485000,  6.01095e-09] 

qo  =  l.lq* 

[51.612000,  1.814500,  -1.712928] 

Final  estimate  for  data  with  5%  noise 

qo  =  0.99(7* 

[46.450800,  1.633500,  3.22126e-09] 

qo  =  1-01(7* 

[47.389200,  1.666500,  -6.6658e-10] 

qo  -  0.95(7* 

[44.574000,  1.567500,  5.56261e-09] 

qo  =  l-05g* 

[49.266000,  1.732500,  -1.061871e-08] 

qo  =  0.9?* 

[42.228000,  1.485000,  6.02759e-09] 

qo  =  l.lq* 

[51.612000,  1.815000,  -1.712024e-08] 

Table  2:  Parameter  estimation  results  for 
q*  =  [k;,  k*,  <]  =  [46.92, 1.65, 1.581139e  -  09] 


We  are  able  to  obtain  reasonable  estimates  for  the  parameters,  especially  At*  and  k*^.  We  have 
difficulty  estimating  the  value  of  /t*,  most  likely  because  the  value  is  small.  In  general,  the 
estimation  error  increases  with  the  error  in  the  initial  guess.  However,  increasing  the  noise 
level  in  the  data  does  not  significantly  effect  the  estimation  accuracy. 

We  note  that  we  obtain  better  estimates  for  the  mean  values  (7o,Co,  and  Tq)  than  for  the 
pressure  coefficients  (Ac*,At^,  and  /c*).  This  is  understandable,  as  the  mean  values  are  more 
influential  in  the  system  dynamics.  They  are  also  more  important  in  identifying  and  char¬ 
acterizing  the  material. 

In  an  electromagnetic  interrogation  parameter  estimation  problem,  an  estimate  is  suffi¬ 
cient  if  it  can  be  used  to  classify  the  material.  We  consider  the  results  for  estimating 
q*  =  [av*,  ac^,  k*\  =  [46.92, 1.65, 1.581139e  —  09]  using  data  with  5%  relative  normal  noise  and 
an  initial  guess  with  -10%  error.  After  solving  the  parameter  estimation  problem,  we  obtain 
the  result  q  =  [42.228000, 1.485000, 6.02759e  —  09].  If  parameter  values  within  the  hypothet- 
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ical  range,  40  <  <  50,  1.4  <  <  1.8,  and  le  -  09  <  /c^  <  9e  -  09  are  characteristic  of 
the  material  under  interrogation,  we  are  successful  in  our  attempt  to  solve  the  estimation 
pioblem.  On  the  other  hand,  if  the  characteristic  material  parameters  fall  within  the  (hy¬ 
pothetical)  range  45  <  k,^  <  47,  l.C  <  <  1.7,  and  le  -  09  <  k,  <  2e  -  09,  we  are  unable 

to  characterize  the  material  with  our  estimates  and  our  attempt  is  unsucce.ssfuL  Ranges  of 
these  parameter  values  for  different  materials  have  not  been  experimentally  determined,  so 
we  have  no  concrete  measure  as  yet  to  asse.ss  our  ability  to  solve  the  problem. 


5  Concluding  remarks 

We  have  presented  theoretical  and  computational  results  for  a  new  class  of  inverse  problems 
arrsing  in  nondestructive  interrogation  of  materials.  Our  focus  is  on  reflections  from  acoustic 
pressure  waves  that  are  moving  through  dielectric  material  targets.  A  detailed  atomic  based 
model  for  acoustic-dependence  of  dielectric  parameters  in  a  Debye  material  was  given  and 
this  was  incorporated  into  a  theoretical  framework  for  both  forward  solutions  and  least 
squares  parameter  estimation. 

Computational  findings  suggest  that  primary  parameters  (relaxation,  static  permittivity, 
etc.)  will  be  readily  identifiable  while  aspects  of  the  nonlinear  dependence  on  pressure  (finst 
order  coefficients)  may  be  ascertained  if  appropriate  data  is  available. 

Our  efforts  on  this  methodology  are  continuing.  An  experimental  device  (similar  to  that 
depicted  in  Figure  7.1  of  [4])  is  currently  under  construction.  Data  from  this  device  will  be 
used  to  test  and  validate  the  methods  developed  in  this  paper. 
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